Quantum Breathers in a Nonlinear Lattice 



W. Z. Wang*, J. Tinka Gammel, A. R. Bishop, 
Theoretical Division, Los Alamos National Laboratory, Los Alamos, New Mexico 87545, USA 



N£> 
On 
On 



On 



> 
oo 

7— I 

m 
O 

On 
c3 



o 
o 



and M. I. Salkola 

Department of Physics & Astronomy, McMaster University, Hamilton, Ontario L8S 4M1, Canada 

(February 1, 2008) 

We study nonlinear phonon excitations in a one-dimensional quantum nonlinear lattice model using 
numerical exact diagonalization. We find that multi-phonon bound states exist as eigenstates which 
are natural counterparts of breather solutions of classical nonlinear systems. In a translationally 
invariant system, these quantum breather states form particle-like bands and are characterized by 
a finite correlation length. The dynamic structure factor has significant intensity for the breather 
states, with a corresponding quenching of the neighboring bands of multi-phonon extended states. 

1995 PACS: 63.20.Ry,63.20.Hp,11.10.Lm,03.65.Ge 



The paradigms of nonlinearity have provided numer- 
ous insights into condensed matter physics ]jj . For exam- 
ple, the dynamics of nonlinear excitations (such as soli- 
tons, polarons, and breathers) are central to understand- 
ing thermodynamic and transport properties of various 
low dimensional materials 0. As an important example 
of intrinsic dynamic nonlinearity, breather excitations — 
spatially localized and time-periodic waves in the form 
of bound states of linear excitations — have been found 
or excluded in various nonlinear models depending on 
the properties of nonintegrability and discreteness HQ- 
In particular, although breathers are rigorously stable 
in integrable partial differential equations (such as the 
(l+l)-dimensional sine-Gordon or nonlinear Schrodinger 
equations), they are unstable in nonintegrable continu- 
ous systems. Recently it has been appreciated, however, 
that they can be stablized, not only in integrable m 
but also in nonintegrable Q cases, by lattice discrete- 
ness and sufficiently strong nonlinearity — in the form 
of dynamic localized excitations, the generalization of un- 
coupled oscillators. For some models, existence regimes 
have been rigorously established ||. We are left with a 
central question for physical problems described by dis- 
crete, nonintegrable quantum systems: do nonlinear solu- 
tions analogous to the classical breather exist and, if so, 
what are their observable signatures? This question has 
a long history in terms of biphonon bound states or res- 
onances H , and has been studied more recently in Bethe 
Ansatz systems |l(| as well as semi-classically |ll]] . Here, 
we focus on numerical exact diagonalization of a simple 
nonintegrable quantum anharmonic chain. We find that 
the classical breather solutions indeed have their quan- 
tum counterparts as eigenstates of the system Hamilto- 
nian with distinctive signatures in appropriate correla- 
tion functions. 

Our quantum nonlinear lattice model is 



H = 7ifi V [- 



$n+l) 2 +tU0n + v 4>i 
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with periodic boundary conditions, and with p n and <p n 
the dimcnsionless canonical lattice momentum and dis- 
placement operators, [4> m ,Pn] = iS mn . The coefficient 
hil sets the scale of energy, whereas n, w, and v are di- 
mensionless parameters. The first two terms describe a 
standard harmonic Debye lattice, while the second two 
terms for w, v>0 describe a single- well </> 4 nonlinearity, 
which corresponds to the Taylor expansion of either a 
local nonlinear interaction due to anharmonicity, or an 
effective lattice interaction arising from, e.g., electron- 
lattice coupling. Other nonlinear Hamiltonians can be 
treated by the same methods described below and can 
be expected to contain varieties of breathers. However, 
the simple form in Eq. (|IJ) isolates the breather with- 
out, e.g., complications of multiply degenerate ground 
states. Note that the class of models represented by Eq. 
(JjJ) does not conserve the total number of phonon quanta 
in contrast with simpler discrete quantum models where 
solitons have been discussed [Q . 

There are two limiting cases for which the physics of 
the model becomes transparent. First, in the absence of 
nonlinearity (v=0), the model is trivially solvable. The 
linear coupling between nearest-neighbor oscillators in 
Eq. (|l|) results in spatially extended Bloch-wave states 
of phonons.The w term acts merely as an adjustable fac- 
tor controlling the dispersion of phonon bands, and will 
be kept small throughout this study. Second, for v (~1) 
3> n, the 4> 4 term plays the dominant role. However, 
we can again qualitatively understand many features of 
the system, including the appearance of breathers. The 
nearest-neighbor interaction governed by n will cause the 
local "Einstein" modes to become hybridized, creating 
narrow continua of phonon states at each local vibra- 
tional energy. The </> 4 term produces a local repulsive 
interaction between two or more quanta if they are lo- 
cated at the same site. This interaction therefore leads 
to bound multi-phonon states, breathers, that are split 
upwards from each multi-phonon continuum pp. These 
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states are the simplest form of breathers we expect to find 
in molecular systems, i.e., whenever local nonlinear vi- 
brations are weakly coupled to neighboring ones. Below, 
we consider systems with 77~u~l, which are more inter- 
esting because of their relevance to many strongly inter- 
acting solid-state systems, such as quasi-one-dimensional 
polymers. In this case, the nearest-neighbor Debye cou- 
pling will compete with the c/> 4 nonlinearity to determine 
the correlation length of eigenstates in the system. The 
main physical effects of the (f> 4 term are: (i) opening finite 
gaps in the many-body excitation spectrum, which offers 
the possibility of bound states within the gaps; and (ii) 
suppression of the effective linear coupling, facilitating 
the dynamic localization of certain eigenstates. 

We study this problem by numerical exact diagonal- 
ization using the Einstein phonon basis |[l3fl , where the 
only physical approximation is the truncation of the in- 
finite Hilbert space. This restricts our numerical simu- 
lations to a parameter region where the nonlinearity is 
not too large. To maintain physical correctness and nu- 
merical reliability, we solve the problem for finite chains 
with moderately large nonlinearity, systematically study- 
ing the effect of truncation and system size. We adopt a 
modified JI3] Davidson method [jl5| to perform the exact 
diagonalization for the large sparse matrix of the system 
Hamiltonian, and then a projective Lanczos method |l6| ] 
to calculate various correlation functions. Eigcnfunction 
convergence has been tested through the Weinstein lower 
bound formula Jl5|, and the symmetries (e.g., the trans- 
lational invariance under periodic boundary conditions) 
were also checked. 

The low-lying eigenspectra for 4- and 8-site nonlinear 
chains with representative parameter sets are shown in 
Fig. ??. Several characteristic features can be noted: 
First, as anticipated above, finite gaps are opened due 
to the effect of the nonlinearity v. Second, states which 
originally belong to the same band in the linear Debye or 
Einstein lattice now split into several bands. While the 
nonlinearity does not allow the ordering of these bands in 
terms of bare phonons of the linear system, one can still 
roughly distinguish different bands by the bare phonon 
number distribution function p a (m) = (cn\S(J2 n b^bn — 
m)\a) and the average number of bare phonons (M a ) = 
Hm m Pa(. m )i where b n = (4> n + ip n )/V2, and \a) is an 
eigenstate of the system with crystal momentum q a and 
energy e a . (We have also used this function to deter- 
mine that our results are not dependent on the trun- 
cation.) In contrast to the Debye lattice, the ground 
state is isolated, with a finite gap to the excited eigen- 
states. The quantity p a (rn) for the first band above the 
ground state shows that it is mainly formed from one 
phonon states (peaked at m=l). There are several dif- 
ferent kinds of higher excited bands of states, includ- 
ing ones with small but finite bandwidths. The emer- 
gence of these isolated bands can be readily interpreted 



as the existence of particle-like states, where the band- 
width is a measure of the inverse particle mass. They 
can also be viewed as bound phonon states. We will ar- 
gue below that these narrow, isolated bands of states are 
quantum breathers. Note that these breather states have 
explicit Bloch-wave translational symmetry with well- 
defined crystal momenta q a associated with their center 
of mass motion. However, they also exhibit a finite lattice 
displacement correlation length (describing the particle- 
like coherence of breathers), as we now discuss. 

To identify breathers in the quantum excitation spec- 
trum, we establish their main characteristics by exam- 
ining key correlation functions. These functions are not 
only important for diagnostic purposes, but are also rel- 
evant to many measurable quantities. One characteristic 
of the classical breather is spatial localization of the enve- 
lope. In the translationally invariant quantum case, the 
dynamic localization of breathers can be measured by the 
corresponding spatial correlation of the displacements at 
sites separated by n lattice spacings: 

f a (n-n\t-t') = (a\4> n ,(t')4> n (t)\a}. (2) 

For the linear system, these correlation functions are 
readily computed and they in general have multiple spa- 
tial and temporal Fourier components reflecting the sys- 
tem's phonon modes propagating with various phase ve- 
locities. However, in a nonlinear system, multi-phonon 
bound states can appear. In Fig. ??, the instantaneous 
spatial correlation function f a (n,0) of a typical state in 
the isolated bands is illustrated. The spatial correlation 
is clearly localized at n = and exhibits the particle-like 
nature of these states, a key feature of a breather. In this 
strongly nonlinear regime, the breathers have a small cor- 
relation length (£ is on the order of a lattice constant). 
The other states form a "continuum" of extended multi- 
phonon states which have significant weight even at large 
n, also illustrated in Fig. 2. 

We investigated temporal coherence in the localized 
state by examining the Fourier transform F a (n, u>) of the 
time-dependent correlation function f a (n,t), 

e/3 + (3) 



calculated using the spectral projection method Jl6| . 
F a (0,oj) for a typical breather state is shown in Fig. ??. 
The dynamics of the localized states exhibit a small num- 
ber of frequencies with significant amplitudes, reflecting 
the anharmonicity of the system. The strongly localized 
breather here can be expected to exhibit multiple inter- 
nal frequencies even classically ^,^,0, especially in the 
strongly nonlinear regime. However, this multi-frequency 
property is not distinctive for breathers - anharmonic 
phonons and multi-phonon states exhibit similar signa- 
tures. 
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These quantum breathers do show experimentally ob- 
servable signatures of their distinctive local structure 
and dynamics. In particular, lattice displacements 
within a few correlation lengths £ are strongly corre- 
lated in the breather and therefore can be analyzed 
by the density-density correlation function, S (r, t) = 
(1/N) I dr'(p(r',0)p(r' + r, i)>, where p(r, t) = *(r - 
nao — £(j) n (t)) is the density operator, ao is the lattice con- 
stant, N is the number of lattice sites, and i—ifi/M^l) 1 / 2 
sets the scale of length for lattice displacements. The 
qualitative behavior of S(r, t) is sensitive to the ratio 
r/£. In principle, these correlations can be probed by 
neutron scattering which directly measures the dy- 
namic structure factor S(q,u>) (the spatial and temporal 
Fourier transform of S(r,t)), given at zero temperature 
by 

S{q, u) = - e « + £ o) f dre~ lqr f dr' 

a J J 

i<0|p(r',0)|a)(a|p(r' + r,0)|0), (4) 

where |0) is the ground state. In the linear case, one 
would see a response in S(q, lo) which traces each multi- 
phonon band dispersion, falling off exponentially in q 
(i.e., the quantum Debye-Waller factor) and decreasing 
algebraically with lo, with regular spacing in lo. The re- 
sult for the nonlinear 4-site system of Fig. ??(a) is il- 
lustrated in Fig. 11. The elastic response of the ground 
state at lo—0 has now gained an exponentially decaying q- 
dependence which is different from the harmonic Debye- 
Waller factor |M . We also note that in this zero temper- 
ature calculation, there is no zero-frequency contribution 
from the breathers flT^l ■ I n addition to the expected large 
low-energy contributions from the ground state and the 
first phonon band, breather excitations are the dominant 
contribution, whereas those contributions from the ex- 
tended multi-phonon "continua" are almost negligible. 
The sum rule, J dcu ooS(q,Lo) = flq 2 £ 2 /2, implies that the 
quenching of the extended-state contribution is consis- 
tent with breathers forming as bound states of nearby 
phonons, and is similar to the transfer mechanism of the 
optical oscillator strength to local modes in the presence 
of electronic bound states. The center of the response is 
shifted to higher q for higher breather bands, consistent 
with the fact that higher lying breathers (binding more 
phonons) are narrower in the classical limit. Another 
clear feature of this nonlinear system revealed by S(q, lo) 
is that the breather bands are irregularly spaced. 

To summarize, we have demonstrated that quantum 
breathers exist as eigenstates of the simple nonlinear lat- 
tice system studied with a renormalized mass correspond- 
ing to a narrow but finite bandwidth. The quantum 
breathers have been identified as robust nonlinear so- 
lutions with their own characteristic spatial localization 
and multi-frequency dynamics. Finally, we predict that 



the dynamic structure factor is strongly enhanced at mo- 
mentum and energy transfers corresponding to breather 
excitations and correspondingly reduced elsewhere. This 
novel result should apply in more general situations and 
provide exciting new possibilities for the identification 
of breathers and other nonlinear excitations by neutron 
scattering. The numerical approach used here is readily 
applied to other nonlinear lattices and to electron-phonon 
coupled models, for which results will be presented else- 
where. 
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Figures in "Quantum Breathers in a Nonlinear Lattice 



Figure 1. Eigenspectrum of the nonlinear system: (a) 4-site chain with the phonon 
basis truncated at 17 phonons per site (i.e., A / l p h=18; 65,536 states), 77=0.5, w=0.2, 
v=2.0; (b) 8-site chain with Ai p h=& (1,679,616 states), 77=0.5, w=0.2, u=0.56. q a labels 
the irreducible representation in the translation group. The dispersion is symmetric in 
q a and only eigenstates with q a > are shown. The solid circles denote breather states. 
Other states are plotted as open diamonds. 
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Figure 2 The spatial correlation function / a (ra, 0) of an 8-site system [Figure 1(b)]: The 
solid line is for the lowest breather state with q = 7r/a , (M 42 ) = 2.83, while the gray 
line is for an extended state with q = 7r/a , (M24) = 3.13. 
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Figure 3. The Fourier transform F a (0,u) g k 3 {<-o) of the temporal correlation function for 
the breather state of Figure 2. The delta function peaks are broadened by a Lorentzian 
for visualization purposes. 




Figure 4. The dynamic structure factor S(q,u) of the system in Figure 1(a) at zero 
temperature and I = a o /10. Note that the equal-height $-peaks of S(q,u) at q = 
are not shown. The delta function peaks along the cu-axis are broadened by a narrow 
Lorentzian. 
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